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In the study of ID nonlinear Hamiltonian lattices, the conserved quantities play an important role 
in determining the actual behavior of heat conduction. Besides the total energy, total momentum 
and total stretch could also be conserved quantities. In microcanonical Hamiltonian dynamics, the 
total energy is always conserved. It was recently argued by Das and Dhar that whenever stretch 
(momentum) is not conserved in a ID model, the momentum (stretch) and energy fields exhibit 
normal diffusion. In this work, we will systematically investigate the stretch diffusions for typical 
ID nonlinear lattices. No clear connection between the conserved quantities and heat conduction 
can be established. The actual situation is more complicated than what Das and Dhar claimed. 

PACS numbers: 05.60.-k,44.10.+i,05.45.-a 


I. INTRODUCTION 

The anomalous heat conduction was first predicted for 
the ID Fermi-Pasta-Ulam /3 (FPU-/3) nonlinear lattices 
by Lepri et al in 1997 [Ij. In this pioneering work, it was 
found numerically that the thermal conductivity k di¬ 
verges with the system size N as k oc with 0 < a < 1 
which breaks the Fourier’s heat conduction law [Ij. Nu¬ 
merical simulations also confirm this anomalous heat con¬ 
duction in diatomic Toda lattice [H, carbon nanotubes 
3 and single polymer chains Q, to name a few. On 
the other hand, the ID nonlinear lattices with external 
on-site potential such as Frenkel-Kontorova (FK) and 
lattices show normal heat conduction Much efforts 

had been devoted to unraveling the physical mechanism 
behind normal and anomalous heat conduction in low di¬ 
mensional systems [sl-l^ . The consensus reached in this 
community is that the momentum conservation and di¬ 
mensionality play the important roles in determining the 
actual heat conduction behavior [43l - l45j| . The mode cou¬ 
pling theory predicts that k oc iV“, k oc ha N and 
K oc const for ID, 2D and 3D momentum conserved sys¬ 
tems, respectively. The numerics in 2D and 3D lattice 
systems were found to be consistent with these predic¬ 
tions [dbl - fsijl . In particular, the prediction of length- 
dependent anomalous heat conduction were also veri¬ 
fied experimentally in ID nanotubes and molecular 
chains and 2D suspended graphene [s^. However, 
there is one exception of ID coupled rotator lattice which 
displays normal heat conduction behavior despite its mo¬ 
mentum conserving nature (hsl - fs^ . 

The traditional numerical methods used to investigate 
the heat conduction problem are Non-Equilibrium Molec- 
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ular Dynamics (NEMD) and equilibrium Green-Kubo 
(GK) method [43,13. ^ novel diffusion method in ther¬ 
mal equilibrium was proposed by Zhao (3 which opens 
a new way to explore the heat transport problem in non¬ 
linear systems [Sa IfiOl ■ The mean square displacement 
of energy diffusion generally follows a power-law time de¬ 
pendence as (Ax'^{t))^ oct^. It has also been rigorously 
proven [blj that this energy diffusion method is equiv¬ 
alent to the Green-Kubo method where the connection 
relation of a = B — 1 firstly proposed from particle diffu¬ 
sion analysis [3 can be derived. 

There are also continued theoretical efforts ranging 
from early mode-coupling theory [1, H, 13 1 renormaliza¬ 
tion group analysis (3 : hydrodynamical theory [2^ 13 ; 
self-consistent mode-coupling theory [3! to recent non¬ 
linear fluctuating hydrodynamical theory [62l - [3 - Al¬ 
though there is still debate about the actual classification 
and divergent exponents of the universal classes, these 
theoretical works have greatly improved our understand¬ 
ing on the nature of the anomalous heat transport in low 
dimensional systems definitely. 

Most recently, it was claimed by Das and Dhar that 
“whenever stretch (momentum) is not conserved in a one¬ 
dimensional model, the momentum (stretch) and energy 
fields exhibit normal diffusion” [3. The ID coupled ro¬ 
tator lattice was taken as the example to support this 
claim. However, after carefully studying some typical ID 
nonlinear lattices with normal heat conduction or energy 
diffusion behaviors, we found no obvious connection be¬ 
tween the stretch and momentum conservation and nor¬ 
mal energy and stretch diffusion can be established. Our 
numerical results indicate that the actual situation might 
be more complicated than what has been claimed. 

This paper will be organized as the following. In Sec¬ 
tion H, we will present the detailed numerical results of 
stretch and energy diffusion for typical ID nonlinear lat¬ 
tices such as coupled rotator, FK, combined (FK-|-(/)^) 
lattices. The conclusions will be summarized in Section 
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FIG. 1: (color online). The 0“^ lattice, (a) The energy 
correlation function pE{i,t) and (c) stretch correlation func¬ 
tion pD{i,t). (b) The MSD of energy and (d) 

of stretch. The energy correlation function pE{i, t) 
follows the Gaussian distribution when correlation time t > 
100 while the stretch correlation function pD{i,t) fails to fol¬ 
low the Gaussian distribution. As a result, the MSD of energy 
depends on time linearly as (Aa;^(t))^ oc t, displaying the 
normal diffusion behavior. In contrast, the MSD (Aa:^(t))^ 
of stretch saturates to a constant value after a short time 
scale. In panel (a), the pE{i,t) is shifted upward with a con¬ 
stant value of 1/(A — 1) to maintain vanishing tail close to the 
boundaries [s^. The energy density is set as E — 0.4 and the 
corresponding temperature is around T « 0.44. The lattice 
size is chosen as N = 801. 




FIG. 2: (color online). The coupled rotator lattice of case (i). 
(a) The energy correlation function pE{i,t) and (b) stretch 
correlation function poii, t). The energy correlation functions 
PE{i,t) follow the Gaussian distributions at correlation times 
t = 100, 300 and 500, implying a normal diffusion behavior. 
However, the stretch correlation functions po(i, t) all collapse 
to the same pattern curve as that at t = 0. The energy density 
is set as E = 0.5 and the corresponding temperature is around 
T « 0.54. The lattice size is chosen as N — 601. 


where AHi{t) = Hi{t) — {Hi) and (•) denotes the en¬ 
semble average or time average equivalently for ergodic 
systems. The stretch distribution function p]j{i,t) can 
also be defined similarly asfb^ 


{Amt)AD^ 
^ {ADo{0)ADom 


(3) 


III. 


II. STRETCH DIFFUSION IN TYPICAL ID 
NONLINEAR LATTICES 


We consider the following Hamiltonian for general ID 
lattices 


H 




+V {qi+1 — Qi) + U{qi) 


( 1 ) 


where qt and pt denote the displacement and momen¬ 
tum for the i-th atom, respectively. The interaction po¬ 
tential V{qi+i — qi) only depends on the displacement 
difference of (qt+i — qi). The existence of on-site po¬ 
tential U{qi) will break the conservation of total mo¬ 
mentum. For simplicity, periodic boundary conditions 
qi = qN+i are applied. The atom index i is assigned as 
— {N— l)/2,..., —1, 0,1,..., {N— l)/2 where an odd num¬ 
ber of lattice sizes N is chosen without loss of generality. 

In order to study the stretch diffusion and energy dif¬ 
fusion behaviors for ID lattices, we adopt the definition 
for the excess energy distribution function pE{i,t) as(^ 


{AH,{t)AHo{0)) 

{AHo{0)AHo{0)) 


( 2 ) 


where the local stretch Di{t) = qi+i{t) — qi{t) and 

AA(i) = A(i)-(A). 

In microcanonical systems with periodic boundary 
conditions, both the total energy H = J^i Hi and total 
stretch D = J^i Di = ~9i) are conserved quanti¬ 

ties. As a result, the excess energy distribution function 
psii, t) and the stretch distribution function poih t) sat¬ 
isfy the sum rules as y).- P F.ii, t) = J2i Poiht) = 0 in 
microcanonical svstems|3^ by noticing that X)- Ai(t) — 
J2i (Ai) = 0 with Aj = Hi, Di. 

One can also define a momentum distribution function 
pp{i,t) [H, [ 5 ^. However, for lattices where momentum 
is not conserved, the summation of pp{i,t) is not time in¬ 
variant as '^i pp{i,t) 7^ 0. In this case, it is meaningless 
to discuss the momentum diffusion. 

In numerical simulations, the fourth order symplectic 
algorithm will be used to integrate the equations of mo¬ 
tions for ID lattices. The time steps At = 0.1 or 0.05 
will be adopted. With this numerical setup, the sum of 
PE{ht) ^nd Poii A) can be maintained within the 
range of the order of 10“® and 10“^^, respectively. The 
energy density E = H/N is the input parameter and 
the temperature T = is a derived quantity as for 
isolated microcanonical systems. 
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FIG. 3: (color online). The stretch correlation functions 
pD{i,t) for coupled rotator lattice in case (ii) and (hi). At 
correlation time t = 10, both pD{i,t) show similar pattern 
for case (ii) and (iii) except that the amplitude in case (iii) 
is much larger than that in case (ii). At t = 30, pD{i,t) still 
maintains a clear pattern for case (iii) while the pattern dis¬ 
appears for case (ii). The parameters are the same as that 
used in Fig. O 


A. lattice 


We first consider the ID cf)'^ lattice with the following 
Hamiltonian 




Pf ^ f ^ 1 4 

Y + - Qi) + 


( 4 ) 


The ID (f>‘^ lattice is a typical nonlinear lattice with on¬ 
site potential which shows normal heat conduction be¬ 
havior 00. The total momentum is not conserved due 
to the existence of external on-site potential. It has been 
verified that the energy diffusion is normal as wellfs^. 

This normal diffusion for energy can be seen from 
Fig. [11(a) and (b). The excess energy distribution func¬ 
tions pE{i,t) collapse to an almost Gaussian distribu- 

i2 

tion pE{i,t) ^ at long enough correla¬ 

tion times. As a result, the Mean Square Displacement 
(MSD) of energy follows a linear time depen¬ 

dence as {Ax^{t))^ 2DEt, asymptotically. Here De 
denotes the diffusion constant for energy. 

However, the stretch distribution poiht) fails to fol¬ 
low the Gaussian distribution, as can be seen in Fig. [T] 
(c). The two humps at correlation time t = 100 spreads 
rapidly over the lattice and disappears at latter correla¬ 
tion times. In Fig. IH the MSD (Aa;^(t))^ of stretch 
is plotted as the function of correlation time t. The 
(Aa;^(t))^ saturates to a constant value after a short 
correlation time scale. It is definitely not the normal dif¬ 
fusion behavior which is predicted by Das and Dhar in 
Ref. [t^. The momentum is not conserved for ID (f>'^ 
lattice, while its stretch diffusion is not normal! 


-O— Case (iii) 
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FIG. 4: (color online). The sum of poii, t) for coupled rotator 
lattice in case (iii). It can be seen that poii, t) decays from 
1 to 0 very quickly which is consistent with the result in Ref. 
M- As a comparison, the sum poii, t) for case (ii) can be 
maintained within the order of in the whole correlation 

time range. This indicates that the stretch is conserved in 
case (ii) but not conserved in case (iii). The parameters are 
the same as that used in Fig. (2] 


B. Coupled rotator lattice 

The ID coupled rotator lattice has the following Hamil¬ 
tonian 


^ = E 


-^ + [I - cos(<7,+i - gO] 


( 5 ) 


Although it conserves the total momentum, it has normal 
heat conduction behavior (ssl . as well as normal en¬ 
ergy diffusion behavior 15811. In particular, its momentum 
diffusion is also normal [58l| . 

The stretch conservation is a tricky issue for coupled 
rotator lattice due to the 27r degeneracy of qi. The dy¬ 
namics of the system is invariant to the arbitrary shift of 
multiple 27r for every qi as qi ^ qi -\- 2mT, where n can 
be any integer number. Depending on how to limit the 
qi, the stretch of coupled rotator lattice can be a con¬ 
served quantity or not. We consider the following three 
limitations for qi or Di. 

(i) No limitations. Nothing is done to the value of qi. 
The qi can take whatever it takes during the evolution of 
the systems. In this situation, the stretch is a conserved 
quantity. The local stretch Di = qi+i — qi can take value 
from negative infinity to positive infinity and the parti¬ 
tion function is not well defined [63, [z3. Although this 
effect will cause problem in theoretical analysis, the dy¬ 
namics of the system will not be affected. The energy 
diffusion is normal as can be seen from Fig. [2] (a). In 
this situation, the stretch correlation functions pD{i,t) 
all collapse to the same pattern as that at f = 0 in Fig. 
[5] (b). This might be due to the unbounded fact of the 
value of qi- 

(ii) The qi is limited within —tt < qi < tt. After each 
time step in numerical simulations, the qi is forced to 
reshifted into this region whenever it jumps out. As a 
result, the local stretch Di lies within —27r < Di < 27r. 
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FIG. 5: (color online). The FK lattice, (a) The energy 
correlation function pE{i,t) and (c) stretch correlation func¬ 
tion (b) The MSD (Aa;^(t))^ of energy and (d) 

of stretch. Both the energy correlation function 
pE{i,t) and the stretch correlation function poiiyt) follow 
the Gaussian distribution after long enough correlation time. 
As a result, both the energy and stretch diffusion is normal 
as oc t asymptotically. In panel (a) and (c), 

the Pe/ oiht) are shifted upward with a constant value of 
1/{N — 1) to maintain vanishing tail close to the boundaries. 
The parameter for on-site coupling strength is set as V = 1. 
The energy density is set as if = 1 and the corresponding 
temperature is around T ~ 0.86. The lattice size is chosen as 
N = 801. 


The stretch is a still conserved quantity and the parti¬ 
tion function is well defined. In numerical simulations, 
the sum "Yhi Poibi) can be maintained within the order 
of 10“^^ for all times which is a signature of conserved 
quantity. The poiijt) displays similar spatial pattern at 
short correlation times as that in case (iii). The only dif¬ 
ference is that the amplitude is much smaller in case (ii), 
as can be seen from Fig. El (a). At larger correlation times 
see Fig. [21(b), the poih t) quickly loses its spatial pattern 
in comparison to that in case (iii). The MSD (Aa;^(t))^ 
of stretch in this case saturates to a constant value after 
a short correlation time (not shown here), similar to that 
of lattice. 

(iii) The local stretch Di is limited within —tt < Di < 
TT. The Qi is not affected during the dynamical evolution. 
However, Di is adjusted appropriately every time when it 
is recorded to generate the correlation function. There¬ 
fore, the partition function is well defined. Only in this 
special situation, the stretch is not a conserved quantity 
as can be seen from Fig. E) 

From the above results and discussions, it is found that 
the conservation of stretch is a very tricky issue. Depends 
on the limitation of qi or Di, the stretch can be adjusted 
to be a conserved or nonconserved quantity. 


FIG. 6: (color online). The combined ¥K+(f>^ lattice, (a) The 
energy correlation function pE{i,t) and (c) stretch correlation 
function pD{i,t). (b) The MSD (Aa;^(t))^ of energy and 
(d) (Aa;^(f))^ of stretch. The overall behavior is similar to 
the (j)^ lattice. Only the energy correlation function pE{i,t) 
follows the Gaussian distribution. As a result, only the energy 
diffusion is normal as (Aa;^(t))^ oc t asymptotically. In panel 
(a), the pE{i,t) is shifted upward with a constant value of 
1/(A — 1) to maintain vanishing tail close to the boundaries. 
The parameter is set as F = 0.5. The energy density is set 
as E — 0.5 and the corresponding temperature is around T « 
0.47. The lattice size is chosen as N = 801. 


C. FK lattice 


We then consider another nonlinear lattice with on-site 
potential, the FK lattice with Hamiltonian as 


^ = E 


r^2 


y + - Qif + y(1 - cos27rgi) 


( 6 ) 


The FK lattice also exhibits normal heat conduction Q 
as well as normal energy diffusion behaviors. In Fig. 
[5] (a) and (b), the excess energy distribution function 
PE(i,t) and the MSD (Ax^(t))^ of energy are plotted. 
The pE{i,t) follows the Gaussian distribution functions 
and the (Ax^(t))^ is linearly proportional to the corre¬ 
lation time as (Aa;^(t))^ oc t, indicating obvious normal 
diffusion behavior for energy. 

In contrast to the cj)'^ lattice, the stretch distribution 
function p£){i,t) approaches the Gaussian distributions 

i2 

as pniijt) ^ -y====e~*^^ at long enough correlation 
times for the FK lattice, as can be seen in Fig. |5| (c). 
The MSD (Ax^(t))^ of stretch do follow the linear time 
dependence as (Aa;^(t))^ 2DDt asymptotically in in 
Fig. |S](d). The stretch diffusion is normal for FK lattice 
despite its stretch conservation nature! 

Unlike the coupled rotator lattice, there is no ambigu¬ 
ous space for the stretch conservation in FK lattice. The 
Qi or Di is not 27r invariant anymore due to the existence 
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of the interaction potential in Eq. 


D. Combined (FK+0^) lattice 


In the end, we consider the combined (FK+^^) lattice 
with Hamiltonian 


^ = E 


y + - (lif + “ cos27rgi) + igf 


(7) 

The combined FK+t/)"^ lattice should also have normal 
heat conduction behavior. This can be verified by ex¬ 
amining the energy diffusion behavior iir Fig. |6] (a) 
and (b). The pE{i,t) follows the Gaussian distribu¬ 
tions and the MSB of energy depends linearly on time 
as (Ax^(t))^ ~ 2DEt asymptotically. 

In Fig. [6](c), the stretch distribution functions poih t) 
are plotted for correlation times t = 100,300 and 500. No 
Gaussian behavior is observed and the scenario is similar 
to that of (j)^ lattice as in Fig. [I](c). Same as cj)'^ lattice, 
the saturates to a constant value after a short 

time scale, as cair be seen in Fig. [5] (d). This is another 
example that the momentum is not conserved, while the 
stretch diffusion is not normal. 


III. CONCLUSIONS 

In conclusion, we have systematically investigated the 
stretch diffusion as well as the energy diffusion for a few 
ID nonlinear lattices with normal heat conduction be¬ 
haviors. For isolated lattices with periodic boundary 


conditions, both the total energy and total stretch are 
conserved quantities. Depends on the existence of on¬ 
site potential, the total momentum can be conserved or 
nonconserved. For (j)'^ and combined (FK-|-^"‘) lattices, 
the total momentum is not conserved while the stretch 
diffusion is not normal. For FK lattice, the total momen¬ 
tum is not conserved and the stretch diffusion is normal. 
For coupled rotator lattice with normal momentum dif¬ 
fusion, the situation is tricky in the sense that its stretch 
conservation depends on the choices of limitation of qi 
or Di. Our numerical results do not support the claim 
that “whenever stretch (momentum) is not conserved in 
a one-dimensional model, the momentum (stretch) and 
energy fields exhibit normal diffusion” proposed by Das 
and Dhar in Ref. 0. However, there is something inter¬ 
esting for the lattices with cosine or bounded potentials. 
It is still an open issue and we hope more efforts will be 
done in this direction in the future. 
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